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Collective transverse flow in heavy-ion collisions at incident energies i?i a b — (1—160)^4 GeV is an- 
alyzed within the model of 3-fluid dynamics (3FD). Simulations are performed with purely hadronic 
equation of state (EoS). At the AGS energies the flow turns out to be sensitive to the stopping 
power of nuclear matter rather than only to the stiffness of the EoS. When the stopping power is 
fixed to reproduce other observables, the flow data favor more and more soft EoS with the incident 
energy rise, which can be associated with "a transition from hadronic to string matter" reported in 
the Hadron- String-Dynamics (HSD) model. Problems, which are met in simultaneous reproduction 
of directed and elliptic flows within the 3FD, suggest that the transverse flow is very sensitive to the 
character of the transverse-momentum noneqmlibrium at the initial stage of collision. Arguments 
in favor of "early-stage" nature of the flow observable are put forward. This suggests that the 
flow (especially the directed one) is determined by early-stage evolution of the collision rather than 
freeze-out stage. 
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I. INTRODUCTION 

The interest to nucleus-nucleus collisions at incident 
energies E\ a ^ — (10-40) A GeV has been recently revived, 
since the highest baryon densities 0,0,13 and highest rel- 
ative strangeness 0, Q at moderate temperatures are ex- 
pected in this energy range. The onset of deconfincment 
is also expected in this domain. In particular, the energy- 
scan SPS program |(J is dedicated to the search for the 
onset of deconfincment in heavy-ion collisions. Moreover, 
a critical end point Q of the QCD phase diagram may 
be accessible in these reactions [H, @. The above ex- 
pectations motivated the project of the new accelerator 
facility FAIR at GSI H3, the heavy -ion program of which 
is precisely dedicated to studying dense baryonic matter 
with the emphasis on the onset of deconfinement and the 
critical end point. The future SPS H3 and RHIC H3 
programs are also devoted to the same problems. 

In Refs. 0, 0] we have introduced 3-fluid dynamical 
(3FD) model which is suitable for simulating heavy-ion 
collisions in the range from AGS to SPS energies, which 
overlaps the energy range of the future FAIR. We have 
started simulations [j| with the purely hadronic equation 
of state (EoS) [l4| in order to see how good the available 
data from AGS and SPS can be understood without in- 
volving the concept of deconfinement. This EoS should 
serve as a reference point for simulations involving more 
sophisticated EoS's including the phase transition. 

With this simple hadronic EoS we succeeded to rea- 
sonably reproduce a great body of experimental data in 
the incident energy range i?i a b — (1-160)A GeV. The list 
includes rapidity distributions, transverse-mass spectra, 



and multiplicities of various hadrons. However, we also 
found out certain problems. In particular, we failed to 
describe the transverse flow at > 40^4 GeV. 

One of the main conclusions of Ref. Q is that in or- 
der to conclude on the relevance of a particular EoS the 
whole available set of data should be analyzed in a wide 
incident energy range with the same fixed parameters of 
the model. Indeed almost for any particular piece of data 
the parameters of the model can be tuned in such a way 
that the model reasonably reproduces the data. At the 
same time the reproduction of the whole set of available 
data with the same model parameters is far from being al- 
ways possible. In Ref. we performed still fragmentary 
analysis of data due to constraints to frames of a single 
paper. Therefore, in the present paper we would like to 
continue our analysis of a particular subset of the data 
concerning collective flow in nuclear collisions, based on 
the set of model parameters fixed in Q. Accordingly to 
our preliminary results 0], this subset is the most prob- 
lematic for the purely hadronic EoS we use. 

The collective transverse flow is formed at the early 
(compression) stage of the collisions and hence reveals the 
early pressure gradients in the evolving nuclear matter. 
The harder EoS is, the stronger pressure is developed. 
Thus, the flow reflects stiffness of the nuclear EoS. The 
collective flow observable has been extensively exploited 
to obtain information on the E oS 11 1 
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In particular, it was expected that the transition to the 
quark-gluon phase results in significant reduction of the 
directed flow (so called "softest point" effect), since 
the pressure in the quark-gluon phase is lower than that 
in the hadronic phase. Calculations within conventional 
(1-fluid) hydrodynamics |30| indicated almost disappear- 
ance of the flow at the onset of the first-order phase tran- 
sition. However, the nuclear finite stopping power essen- 
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tially weakens this effect |3J, |42j . In the present paper we 
are going to discuss the transverse flow in terms of stiff- 
ness of the nuclear EoS and the stopping power achieved 
in the collision. 

At present there are quite systematic collective-flow 
data, both on directed and elliptic transverse flows, at 
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energies. We confine our present discussion to the flow 
of non-strange probes, i.e. nuclcons and pions. Recent 
analysis of this set of data was done in Ref. 0] with em- 
phasis on momentum dependence of nuclear forces. This 
was done within the JAM model [64| combined with a co- 
variant prescription of mean fields (RQMD/S) [H3 with 
the emphasis on the momentum dependence of nuclear 
forces. We would like to analyze the same data but within 
the hydrodynamic approach. 



II. 3FD MODEL 

A direct way to address thermodynamic properties of 
the matter produced in these reactions consists in appli- 
cation of hydrodynamic simulations to nuclear collisions. 
However, finite nuclear stopping power, revealing itself 
at high incident energies, makes the collision dynamics 
non-equilibrium and prevents us from application of con- 
ventional hydrodynamics especially at the initial stage of 
the reaction. Since the resulting non-equilibrium is quite 
strong, introduction of viscosity and thermal conductiv- 
ity does not help to overcome this difficulty, because by 
definition they are suitable for weak non-equilibrium. A 
possible way out is taking advantage of a multi-fluid ap- 
proximation to heavy- ion collisions. 

Unlike the conventional hydrodynamics, where local 
instantaneous stopping of projectile and target matter 
is assumed, a specific feature of the dynamic 3-fluid de- 
scription is a finite stopping power resulting in a counter- 
streaming regime of leading baryon-rich matter. Exper- 
imental rapidity distributions in nucleus-nucleus colli- 
sions support this counter-streaming behavior, which can 
be observed for incident energies between few and 200A 
GeV. The basic idea of a 3-fluid approximation to heavy- 
ion collisions 0, 0] is that at each space-time point 
x = {t, x) the generally nonequilibrium distribution func- 
tion of baryon-rich matter, can be represented as a sum of 
two distinct contributions, fhar.( x iP) — fp(. x iP)+ft(x,p), 
initially associated with constituent nucleons of the pro- 
jectile (p) and target (t) nuclei. In addition, newly pro- 
duced particles, populating the mid-rapidity region, are 
associated with a fireball (f) fluid described by the distri- 
bution function ff(x,p). It is assumed that constituents 
within each distribution are locally equilibrated, both 
thermodynamically and chemically. This assumption, 
justifying the term "fluids" , relies on the fact that intra- 
fluid collisions are much more efficient in driving a system 
to equilibrium than inter-fluid interactions. Therefore, 
the 3-fluid approximation is a minimal way to simulate 
the finite stopping power at high incident energies. Note 



that both the baryon-rich and fireball fluids may con- 
sist of any type of hadrons and/or partons (quarks and 
gluons), rather than only nucleons and pions. 

Our 3FD model 0, El is a straightforward extension 
of the 2-fluid model with radiation of direct pions [63.1671 
HI and (2+l)-fluid model [gjllzjj. We extend the above 
models in such a way that the created baryon-free fluid 
(which we call a "fireball" fluid, following the Frankfurt 
group) is treated on equal footing with the baryon-rich 
ones. We allow a certain formation time for the fireball 
fluid, during which the matter of the fluid propagates 
without interactions. We assume that the fireball matter 
gets quickly thermalized after its formation. The latter 
approximation is an enforced one, since we deal with the 
hydrodynamics rather than with kinetics. 
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FIG. 1: (Color online) Baryon-density dependence of the 
pressure at zero temperature for hard (K = 380 MeV), in- 
termediate (K = 210 MeV), soft (K = 130 MeV), extrasoft 
(K — 100 MeV), and gas EoS's. Shaded region is the con- 
straint derived from experimental data in Ref. 

The main unknowns of the present simulations can be 
briefly summarized as follows: an EoS and "cross sec- 
tions". Our goal is to find an EoS which in the best 
way reproduces the largest body of available observables. 
The "cross sections" arc equally important. They de- 
termine friction forces between fluids and hence the nu- 
clear stopping power. In principle, friction forces are EoS 
dependent, because medium modifications, providing a 
nontrivial EoS, also modify cross sections, and should be 
externally supplied together with the EoS. However, at 
present we have at our disposal only a rough estimate 
of the friction forces [fjj- Therefore, we have to fit the 
friction forces to the stopping power observed in proton 
rapidity distributions. 

We have started our simulations |3j with a simple, 
purely hadronic EoS 01 which involves only a density 
dependent mean field providing saturation of cold nu- 
clear matter at normal nuclear density no = 0.15 fm -3 
and with the proper binding energy -16 MeV, cf. Ap- 
pendix El This EoS is a natural reference point for any 
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other more elaborate EoS. 

Fig. ^ displays the density dependence of the pressure 
of four versions of this hadronic EoS at zero temperature: 
hard EoS with incompressibility K = 380 MeV, interme- 
diate EoS with K = 210 MeV, soft EoS with K = 130 
MeV, and extrasoft EoS with K = 100 MeV. The pres- 
sure of the nucleon gas is also presented, indicating the 
softest EoS within our parametrization. As seen from 
Fig. ^ the intermediate EoS, with which all the calcula- 
tions of Ref. were done, and the soft EoS, are within 
the constraint given by Danielewicz et al. [44| based on 
the analysis of flow of nuclear matter at the AGS incident 
energies. 




Vs NN [GeV] 

FIG. 2: (Color online) Incident energy dependence of midra- 
pidity yield of various hadrons produced in central nucleus- 
nucleus collisions. At the AGS energies the 3FD calculations 
with intermediate EoS (K — 210 MeV) were done for Au+Au 
collisions at impact parameter b — 2 fm, at SPS energies — 
for Pb+Pb at b — 2.5 fm. The compilation of experimental 
data is taken from Ref. [TsT ]. 

Strictly speaking, the constraint re gion of Ref. |S| is 
valid the class of EoS's considered in |44|, since the flow 
essentially depends also on the temperature dependence 
of the pressure, which is not displayed in Fig. More- 
over, since the incompressibility characterizes the EoS 
only in the vicinity of the bound state, different EoS's 
with the same incompressibility may differ at higher den- 
sities even at zero temperature. This is the reason why 
the incompressibility of our intermediate EoS (K = 210 
MeV) differs from that of the intermediate EoS of Ref. 

(K = 300 MeV), while these two EoS's are very sim- 
ilar at densities higher than 2hq. The same concerns soft 
EoS's: K = 130 MeV in our case, and K = 210 MeV 
in Ref. 0. Of course, values K = 100 and 130 MeV 
are unrcalistically small from point of view of low-energy 
nuclear physics. However, in simulating heavy- ion col- 
lisions we predominantly deal with densities n > Itiq. 
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FIG. 3: (Color online) The same as in Fig. |5|but zoomed and 
for three EoS's displayed in Fig. Q 

Therefore, we should estimate the degree of reasonable- 
ness of a EoS precisely in this density region, considering 
the incompressibility only as label for the EoS. 

A comprehensive analysis of the nuclear EoS [7^ . 
based on astrophysical data, indicated that the prefer- 
able region of Fig. should be slightly shifted to higher 
pressures. Thus, the hard version of our EoS also falls 
into the region of EoS's reasonable from the astrophysical 
point of view. At the same time, our soft and extrasoft 
EoS's fall out of this astrophysically acceptable region. 

The 3FD model with the intermediate EoS turned 
out to be able to reasonably reproduce a great body 
of experimental data Q in a wide energy range from 
AGS to SPS. Fig. |21 illustrates (and extends) the re- 
sults of Ref. . The experimental data were com- 
piled in Ref. |73| based on reported experimental results 
111 El EI Iz31zl El IH M, HI- These data slightly 
differ in degree of centrality, at which they were taken. 
We however performed our calculations for Au+Au col- 
lisions at fixed impact parameter b = 2 fm for the AGS 
energies, and for Pb+Pb at b = 2.5 fm for the SPS en- 
ergies. Slight irregularity of calculated curves between 
y/SNN = 4.85 GeV and y / s]v7V = 6.27 GeV (i.e. between 
-Blab = 10 A GeV and E lah = 20 A GeV) is related to 
transition to different colliding nuclei at different impact 
parameter. Poor description of K + and A yields at low 
incident energies is not surprising, since our EoS is based 
on the grand canonical ensemble. Really surprising is the 
fact that K~ and A yields are still well reproduced at low 
energies. 

In the present paper we also do calculations with hard, 
soft and estrasoft EoS's (cf. Fig. ^) in order to study ef- 
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fects of the EoS stiffness. Changes, introduced by these 
hard and soft EoS's, are presented in Fig. |21 As seen, 
the intermediate EoS is indeed preferable in overall re- 
production of the data. 

To study effect of nuclear stopping power on the trans- 
verse flow, we also do calculations within the conventional 
1-fluid dynamics (1FD), for which we have a separate 
code. 



III. DIRECTED FLOW 

The conventional transverse-momentum flow of an a 
species is defined as [&| 



J d 2 p T p x E dN a /d 3 p 



d 2 p T E dN a /d 3 p 



(1) 



where p x is the transverse momentum of a particle in the 
reaction plane, E dN a /d 3 p is the momentum distribution 
of a-hadrons, which takes due account of feed-down from 
resonance decays, and integration runs over the trans- 
verse momentum pt- 
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FIG. 4: (Color online) Directed flow of nucleons as a function 
of rapidity for mid-central Au+Au collisions at -Ei a b = 1.15/1 
GeV. The 3FD results at b = 6 fm are presented for three 
EoS's. The data are from Ref. |2l||. 

In Figs. the comparison with Bevalac and AGS 
data on the directed flow of protons [U EE an d "nu- 
cleons" |2l], is presented. Under "nucleon data" we 
mean the combined flow of identified protons and light 
nuclear fragments [2lL IHfij. In fact, we compute (p x ) of 
so-called primordial nucleons, which later may coalesce, 
forming light fragments. In view of this, the comparison 
with the nucleon data is preferable for our model, since 
it avoids additional assumptions of coalescence. In the 
case of -Eiab = 10. 5A GeV (cf. Fig. 0), agreement with 
nucleon data indeed seems to be better than that with 



identificd-proton data. Since experimentally the side- 
ward flow is observed for mid-central collisions with an 
average impact parameter of b = 6 fm, we restrict our 
computation to this single impact parameter here. 

It is seen that the observed directed flow prefers more 
and more soft EoS with the incident energy rise. While 
at £\ a b < 2A GeV the intermediate EoS is preferable, 
even the soft EoS looks too hard at £q a b > 8A GeV. 
At -Ei a b = 8A GeV the extrasoft EoS is certainly prefer- 
able, while at E\ a t, = 10. 5/1 GeV it is difficult to choose 
between soft and extrasoft EoS's. However, in Rcfs. 
[43L liH l48l | , a rather weak sensitivity of the flow to the 
stiffness of the EoS was reported. An apparent reason 
for this disagreement is that in above models the EoS 
becomes effectively softer with incident energy rise. At 
higher incident energies, this softeningoccurs because of 
gradual (governed by form-factors) |43| | or abrupt |48l | 
switching off mean-field interactions, or because more 
and more initial baryon-baryon collisions at the early col- 
lision stage end up in strings which are not affected by 
the nuclear mean field j^. Since the mean-field switch- 
ing off is accompanied by dominating string excitations 
also in models |2J > it is possible to summarize this 
mechanism as "a transition from hadronic to string mat- 
ter" , following Ref. ^3] • This early-stage "string matter" 
results in effective softening of the EoS. 

The flow magnitude is conventionally defined as the 
slope of directed flow at midrapidity 
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The excitation function of the proton flow magnitude is 
shown in Fig. This figure summarizes observations 
made above. In fact, the point at E\ a ^ = 1.15^4 GeV 
is somewhat misleading. While the differential directed 
flow is well reproduced at this energy, cf. Fig. the slope 
is not, because of asymmetry of the data with respect to 
midrapidity. Note that for the symmetric system like 
Au+Au the flow should be antisymmetric with respect 
to midrapidity. Additionally Fig. [7] demonstrates that 
the nuclear stopping power also essentially affects the di- 
rected flow. The conventional 1-fluid simulations (infinite 
stopping power) result in essentially stronger flow than 
it is observed. 

At high energies the azimuthal asymmetry is usually 
characterized by the first coefficients of Fourier expan- 
sion of the azimuthal-anglc dependence of the single- 
particle distribution function, i.e. by the directed flow 
vi = (cos4>) |8l : 



v[ a) (y) 



d 2 PT (Px/Pt) E dN a /d 3 p 



d 2 p T E dN a /d 3 p 



(3) 



Examples of V\ flow for protons and pions at E\ a ^ = 
40A and 158/1 GeV are presented in Fig. [S] for mid- 
ccntral collisions. The 3FD simulations arc performed at 
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FIG. 5: (Color online) Directed flow of protons as a function of rapidity for mid-central Au+Au collisions at E\ a -^ = 2A, 4A, 6A 
and 8A GeV. The 3FD results at b = 6 fm are presented for three EoS's. The data are from Ref. |5S| . 
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FIG. 6: (Color online) Directed flow of nucleons and protons 
as a function of rapidity for mid-central Au+Au collisions at 
^lab = W.5A GeV. The 3FD calculations with three EoS's 
are presented for impact parameter 6=6 fm. Data are from 
Ref. [55|l . Full symbols correspond to measured data, open 
ones are those reflected with respect to the midrapidity point. 
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FIG. 7: (Color online) Proton flow magnitude as a function 
of beam energy for mid-central Au+Au collisions. The 3FD 
calculations are presented for and 6=6 fm. Compilation of 
the data is from Ref. 1581 . 



fixed impact parameter as b = 5.6 fm. Several sets of 
data, which noticeably differ from each other, are shown 
in these figures. The "standard" method [Si IM. H| 
for evaluating the flow coefficients requires an event-by- 



event estimate of the reaction plane with which outgoing 
hadrons correlate. However, this method does not dis- 
criminate other sources of correlations, like those due to 
global momentum conservation, resonance decays, etc. A 
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FIG. 8: (Color online) Directed flow of protons (upper panels) and charged pions (lower panels) in mid-central Pb+Pb collisions 
a t -^lab = 158A (left panels) and 40^1 (right panels) GeV as a function of rapidity. In the upper block of panels (a) the 3FD 
calculations at b — 5.6 fm are presented for intermediate, soft and extrasoft EoS's; at the lower block (b) — the gas EoS and 
EoS with K = 90 MeV. Experimental data Ref. obtained by two different methods are displayed: by the standard method 
(v(st)) and by the method of n-particle correlations (v(n)). Full symbols correspond to measured data, while open symbols 
are those reflected with respect to the midrapidity. Updated data of the NA49 Collaboration loot (v(st) — 98) with acceptance 
0.05 < pr < 0.35 GeV/c for pions and 0.6 < pr < 2.0 Gev/c for protons, are also shown. 
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new method |86( of n-particle correlations allows to get 
rid of these non-flow correlations in extracting V\ and v 2 
from genuine azimuthal correlations. 

At these energies the proton vi flow, calculated with 
intermediate and soft EoS's, significantly differs from the 
data. The experimental data favor extrasoft EoS, which 
became already preferable at £i a b = 8A GeV, cf. Fig. 
Thus the EoS softening with the incident energy rise 
stops near the top AGS energy, at higher energies the 
same extrasoft EoS remains preferable. This fits the ob- 
servation in Ref. 01 that the "transition from hadronic 
to string matter" , i.e. early-stage softening of the EoS, 
practically saturates at the top AGS energy. Apparently, 
the early-stage softening of the EoS at SPS energies is 
the reason of the success of microscopic RQMD |47j . 
UrQMD H3, JAM E3 and HSD models. These 
models reveal quite soft dynamics at the early stage, since 
mean fields turn out to be essentially reduced in this en- 
ergy range. 

The lower block of panels (b) in Fig. [S] demonstrates 
what happens if we further soften the EoS. Slight soften- 
ing of the extrasoft EoS, i.e. taking K = 90 McV, results 
in better reproduction of pion v\ however on the expense 
of proton flow. The softest EoS, i.e. the gas EoS (cf. 
Fig- Q)i already again contradicts the data, since it pro- 
duces a wiggle in proton v±, which is not experimentally 
observed. 

The calculated pion directed flow closely follows the 
pattern of the proton one, while the pion v\ data re- 
veal anticorrclation with proton flow. In fact, our pion 
directed flow correlates with proton one even at lower, 
AGS energies, instead of anticorrclating. A probable rea- 
son for this poor reproduction of pion v\ consists in dis- 
regarding the fact that a part of frozen-out particles is 
"shadowed" by still hydrodynamically evolving matter. 
This mechanism of proton-pion anticorrelation was dis- 
cussed in Refs. |g3> In terms of hydrodynamics, this 
shadowing means that frozen-out particles cannot freely 
propagate through the region still occupied by the hy- 
drodynamically evolved matter but rather get reabsorbed 
into the hydrodynamic phase. Apparently, the baryon di- 
rected flow is less affected by this shadowing. The reason 
is that the baryon directed flow reveals the collective flow 
of matter, since this matter is mainly built of baryons as 
the most abundant and heavy component of the system. 
This collective flow is mainly formed at the early stage 
of the reaction. Baryon rescatterings within this earlier- 
formed collective flow at later stages do not essentially 
alter the collective transverse momentum of the matter. 
At the same time, the pions are much stronger affected 
by this shadowing, since they are screened by the pre- 
dominantly baryonic matter, where pions may rescatter 
or even be absorbed. This can drastically change the 
pion v\. 



IV. ELLIPTIC FLOW 

The elliptic flow, v 2 — (cos 20) |g4|, is the second co- 
efficient of Fourier expansion of the azimuthal-angle de- 
pendence of the single-particle distribution function 



d 2 p T [(p 2 x -p 2 y )/p 2 T ]E dN a /d 3 p 

vi a) (y) = 1 ? ■ (4) 

d 2 p T E dN a /d 3 p 
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FIG. 9: (Color online) Midrapidity elliptic flow as a function 
of impact parameter in Au+Au collisions at E\ a ^ = 2A, 4A 
and 6A GeV. Experimental data for pr > are from Ref. 

Impact parameter dependence of the proton elliptic 
flow for Au+Au collisions at energies E\ a b = 2A, 4A and 
6A GeV without any cut on the transverse momentum 
(pt > 0) is presented in Fig. |3J We see that the soft 
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FIG. 10: (Color online) Proton rapidity spectra at AGS en- 
ergies for various impact parameters. Solid lines correspond 
to 3FD calculations with intermediate EoS. For comparison 
the 1FD result with the same EoS is shown by the dotted- 
dashed line. Experimental points are taken from at 2A 
and 4A GeV, and Ftl at 6A GeV. The percentage indicates 
the fraction of the total reaction cross section, corresponding 
to experimentally selected events. 



and intermediate EoS's, which are the best in reproduc- 
ing the directed flow, are completely off the experimental 
data. At the same time, the hard EoS reasonably repro- 
duces V2- This observation is similar to that done in Rcfs. 
In Refs. HI E [53 it was found that a 
proper momentum dependence in the nuclear mean field 
is of prime importance for the simultaneous reproduction 
of the directed and elliptic flows. 

The fact that this momentum dependence is vitally im- 
portant apparently again implies that the initial stage of 
nuclear collision is essentially noncquilibrium and that 
the flow is sensitive to this early noncquilibrium in 



transverse-momentum distribution. Indeed, even if a mo- 
mentum dependent forces are present in the model un- 
derlying a EoS, in the EoS this momentum dependence is 
integrated out with the equilibrium momentum distribu- 
tion. Therefore, the in- and out-of-reaction-plane trans- 
verse pressures, which are responsible for V\ and V2 flows, 
are tightly interrelated. This prevents the 1FD model 
from simultaneous reproduction of v\ and V2 ■ In the 3FD 
model these two components of pressure are less tightly 
interrelated but still do not possess the same degree of 
freedom as in genuine noncquilibrium, since a superposi- 
tion of three equilibrium distributions certainly confines 
the class of possible noncquilibrium momentum configu- 
rations. The 3-fluid approximation was designed to pri- 
marily simulate the transverse-longitudinal anisotropy of 
momentum distribution. It was not specially tuned to de- 
scribe anisotropy within the transverse direction, which 
is of prime importance for v^.. Apparently, the genuine 
momentum nonequilibrium at the initial stage together 
with proper momentum dependence in the nuclear mean 
field is required for simultaneous reproduction of the di- 
rected and elliptic flows. 

In view of above said, our intermediate EoS, which is 
the best in overall reproduction of a large body of other 
data, can be considered as a kind of compromise allowing 
to simultaneous reproduce orders of magnitude of the 
directed and elliptic flows within the assumptions of the 
3FD model. 
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FIG. 11: (Color online) Pion elliptic flow as a function of 
rapidity in Au+Au collisions at = 11. 5 A GeV. Experi- 
mental data for pr > are from Ref. [Hrl . 



Results of 1FD simulations are also presented in Fig. 
In fact, at energies E lah = 2A,AA and 6A GeV the 
degree of stopping is rather high already in the 3FD. As 
it is seen from Fig. I1UI proton rapidity distributions in 
1FD do not differ v ery much from those in 3FD and from 
experimental data [7a, |7y| . The same concerns other ob- 
scrvables except flow. The flow reveals much harder pat- 
tern within the 1FD at the same intermediate EoS, cf. 
also Fig. [3 and effectively looks like that with the hard 



9 



0.1 
0.05 

*■ 



-0.05 
0.05 



-0.05 



1 1 1 1 1 1 1 1 1 1 T" 

Pb+Pb (158A GeV) 




- P 



H ' — h 



A 


v 2 (st) 


♦ 


v 2 (2) 


■ 


v 2 (4) 


• 


v 2 (st)-98 


1 i 1 i 



1 1 1 1 1 1 1 1 1 1 T" 



Pb+Pb (40A GeV) 



H — ' — 1 



3FD, b=5.6 fm 
— — gas EoS 

ex.soft EoS 

soft EoS 

int. EoS 




_i i i i i i i_ 



A 


v 2 (st) 


♦ 


v 2 (2) 


■ 


v 2 (4) 


• 


v 2 (st)-98 


I 




-3-2-10 1 2 3-3-2-10 1 2 3 

yc.m. yc.m. 

FIG. 12: (Color online) The same as in Fig. [S]but for the elliptic flow. 
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EoS in 3FD. 

Available data on the pion elliptic flow at the AGS 
energies, cf. Fig. ^] are not very spectacular. They 
allow only to estimate the order of magnitude of the vi 
value. As seen, here the difference between hard and 
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FIG. 14: (Color online) The same as in Fig. I13l but for pions. 



intermediate EoS's is not already that large, and the 1FD 
result is in the best agreement with the data. 

Rapidity dependence of the V2 flow for protons and 
pions at i?i a b = 40^4 and 158A GeV are presented in Fig. 
IT21 Here we see that stiffness of the EoS only weakly 
affects the result. Even the difference between extremes, 
i.e. the hard EoS and gas EoS, is quite moderate. The 
same feature was observed in the JAM simulations |45j . 
but the reason there was different. In Ref. ji^i string 
excitations dominate particle production at early times. 
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These strings are not affected by the nuclear mean field 
by the assumption of the JAM model. Therefore, the 
mean-field effects arc small at the SPS energies. In our 
case the mean-field effects (i.e. stiffness of the EoS) are 
present to the full extent, and nevertheless their effect is 
small. 

The pion elliptic flow at 158 A GcV is in surprisingly 
good agreement with experimental data both in mag- 
nitude and shape. However, in spite of the weak de- 
pendence on the EoS stiffness, we still essentially over- 
estimate the proton v-i at both incident energies, 40A 
and 158A GcV. Moreover, the better we reproduce v\ 
(with extrasoft EoS), the worse result we obtain for V2- 
This again implies that the nonequilibrium transverse- 
momentum anisotropy at the initial stage of nuclear 
collision (now already without serious interplay with 
momentum-dependent forces) is very important here. An 
additional argument in favor of this conclusion is that 
cascade simulations j^E E reproduce the proton V2 
quite successfully at 158 A GeV, although still overesti- 
mating at 40^4 GeV. 

The case of energy of 40A GeV deserves special at- 
tention. All the calculations done up to now |45l . Ifilj . 
in JAM, UrQMD and HSD, do not reproduce the "col- 
lapse" of i>2 at 40j4 GeV, which can be interpreted as 
the onset of first-order phase transition (5j|. As for the 
3FD model, the disagreement with the data is the largest 
at 40A GeV, thus also indicating that something special 
happens in the system at this energy. 

The earlier pion V2 data of the NA49 collaboration, 
taken at smaller acceptance at the top SPS energy [6fl |. 
have been already analyzed within hydrodynamic ap- 
proaches. In the expansion model j8^ the Bjorken scaling 
solution |90l | was assumed for longitudinal evolution and 
2D hydro was solved numerically for transverse one. In 
this way, the elliptic flow could be estimated only at the 
midrapidity point. The full 3D expansion model with 
postulated initial conditions was applied to the meson 
elliptic flow by Hirano |46i |. In a qualitative agreement 
with Ref . |89j , it was found that p- meson decays result in 
almost vanishing azimuthal anisotropy of pions near the 
midrapidity. Note that decays of all relevant resonances 
are taken into account in our model. 

Figs. 1131 and ITT1 summarize our results on V2- Be- 
sides, they show that i>2 becomes low sensitive not only 
to the stiffness of the EoS but also to the nuclear stop- 
ping power at high incident energies. Indeed, the 1FD 
results become very similar to that of the 3FD model 1 , 
especially for protons. However, all these results overes- 
timate the proton x>2 at high incident energies. As seen, 
the 3FD with hard EoS and 1FD properly reproduce the 
change of sign of proton V2- Unfortunately, data on pion 
V2 are practically absent at AGS energies. 



1 Note however that the 1FD results are completely unreasonable 
for other observables at high incident energies. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we presented the analysis the collective 
transverse flow of nuclear matter within the 3FD model 
0, El- This is not a dedicated study of flow, aiming 
to reproduce experimental data by fitting model param- 
eters. On the contrary, we have used model parameters 
of our previous paper [3j, which were fixed to reproduce 
a great body of experimental data in the incident energy 
range E\ a b ~ (1-160)^4 GcV. We have varied these pa- 
rameters only in order to study sensitivity of the flow to 
stiffness of the EoS and to the nuclear stopping power. 

This paper, as well as the previous one 0, is based 
on a simple purely hadronic EoS [l4|. We do not expect 
that this EoS perfectly reproduces all the data (and, in 
fact, it does not). However, we believe that some use- 
ful conclusions can be drawn from success or failures of 
simulations based on this simple EoS. 

At the AGS energies the flow turns out to be sensitive 
to the stopping power of nuclear matter rather than only 
to the stiffness of the EoS. The stronger stopping power 
is, the harder dynamics is, i.e. it looks as if a harder EoS 
is employed. This was demonstrated by comparing the 
3FD and 1FD calculations. 

When the stopping power is fixed to reproduce other 
observables, the flow data favor more and more soft EoS 
with the incident energy rise. This EoS softening occurs 
in the AGS energy range, at higher (SPS) energies the 
same extrasoft EoS remains preferable. This observa- 
tion indirectly agrees with that obtained in microscopic 
RQMD [H, UrQMD 0, JAM and HSD jUl mod- 
els. In these models an effective softening of the EoS 
occurs at the early stage of the collision, because of grad- 
ual (governed by form-factors) [U or abrupt switch- 
ing off mean-field interactions, or since mean fields turn 
out to be essentially reduced because of string formation 
miE3- With the incident energy rise, string excitations 
and mean-field switching off result in progressing soften- 
ing of the EoS. With this mechanism of effective soften- 
ing these models efficiently reproduce the directed flow in 
the whole incident-energy range with the same mean-field 
interaction. Since the mean-field switching off is accom- 
panied by dominating string excitations also in models 
|48l l5lj , it is possible to summarize this mechanism as 
"a transition from hadronic to string matter" , following 
Ref. Moreover, as it was stated in Ref. ^(Jj the 

"transition from hadronic to string matter" practically 
saturates at the top AGS energy. This also is in agree- 
ment with our observation that no extra softening of the 
EoS in required at SPS energies, as compared with the 
top AGS energy. 

We have also found that it is impossible to simultane- 
ously reproduce the directed and elliptic flow with the 
same EoS. The directed flow requires a softer EoS, while 
the elliptic flow demands for a harder EoS. This obser- 
vation is similar to that done in Rcfs. [24l lil ]. In Rcfs. 
|43l l45| it was found that a proper momentum depen- 
dence in the nuclear mean field is of prime importance for 
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FIG. 15: (Color online) Directed (four upper panels) and elliptic (four lower panels) flow of protons and charged pions in 
mid-central Pb+Pb collisions at £) a b = 158^4 GeV (left panels) and 40A GeV (right panels) as a function of rapidity. 3FD 
calculations with intermediate EoS at b = 5.6 fm are presented for three different freeze-out criteria: the freeze-out energy 
density £f rz = 0.2 GeV/fm 3 (late freeze-out), £f rz = 0.4 GeV/fm 3 (our conventional freeze-out), and £f rz = 0.6 GeV/fm 3 (early 
freeze-out). 



the simultaneous reproduction of the directed and elliptic 
flows. The fact that this momentum dependence is vi- 
tally important apparently suggests that the initial stage 
of nuclear collision (i.e. the formation of so-called "initial 
fireball") is essentially nonequilibrium and that the flow 
is sensitive to this early nonequilibrium in transverse- 
momentum distribution. 

At SPS energies the directed-flow data favor the same 
extrasoft EoS which was the best beginning from the 
incident energy of 8A GeV. This EoS softening can be 
associated with either the "transition from hadronic to 
string matter" j which has been mentioned already, or 
could be a signal of deconfinement transition. Within the 
3FD simulations the elliptic flow becomes low sensitive 
to the stiffness of the EoS and even to the stopping power 
at these energies. Lack of the proper description of the 



nonequilibrium transverse-momentum anisotropy at the 
initial stage of nuclear collision prevents the 3FD from 
good reproduction of the proton elliptic flow at SPS en- 
ergies. Disagreement of the 3FD with the data is largest 
at 4QA GeV incident energy, which may indicate that 
something special happens in the system at this energy. 

The above discussion suggests that the flow is an 
"early-stage" observable, i.e. that determined by early- 
stage evolution of the collision. Our calculations with 
different freeze-out criteria (different freeze-out energy 
densities £f rz ), displayed in Fig. 1151 confirm this conjec- 
ture. From the point of view of freeze-out, the difference 
between £f rz = 0.2 and 0.6 GeV/fm 3 is huge — they cor- 
respond to completely different patterns of the freeze-out. 
Nevertheless, the directed flow turns out to be fairly in- 
sensitive to this difference. Slight sensitivity is observed 
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only in spectator regions. The elliptic flow is a more sub- 
tle quantity. It is a measure of difference of the in-plane 
and out-of-plane flow. Taking into account that our cal- 
culated vi are of the order of 5% and that sensitivity 
of v\ with respect to above £f rz variations is 1%, the 
value of 20% for the variation of «2 under above change 
of Sf IZ looks quite natural. In fact, this low sensitivity of 
the flow to the freeze-out stage is not surprising. The v\ 
flow is a measure of collective momentum accumulated 
by matter during the expansion stage. The driving force 
of this collective momentum is the pressure gradient cre- 
ated at the early compression stage of the collision. If 
the freeze-out occurs not too early (e.g., not right after 
the compression stage), this pressure gradient has enough 
time to accelerate the matter, and hence late-stage evolu- 
tion does not noticeably change the earlier-accumulated 
collective momentum. 

Thus, problems in reproduction of flow data appar- 
ently result from not quite accurate description of early- 
stage transverse-momentum nonequilibrium within the 
3FD rather than are caused by our intermediate hadronic 
EoS. Therefore, keeping in mind that with our intermedi- 
ate hadronic EoS we have succeeded to reasonably repro- 
duce a large body of experimental data in the incident 
energy range Ei a b ~ (1-160)^4 GeV 0], we can still con- 
clude that this EoS is certainly good (however, not per- 
fect) in this energy range. However, since the required 
EoS softening can be still associated with a deconfmc- 
ment transition, other, more sophisticated EoS's, includ- 
ing phase transition to the quark-gluon phase, should be 
tested within 3FD simulations. The present calculations 
provide a natural benchmark for the future analysis. 
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APPENDIX A: HADRONIC EOS 

The EoS, used in present simulations, was originally 
proposed in Ref. |l4j . The energy density and pressure 



are constructed as follows: 

e(n B ,T) = e gas {n Bl T) + W{n B ), (Al) 

P(n B ,T) = P Eaa {n Bl T)+n B dW } nB) - W(n B ),(A2) 

dn B 

where e gas (n B ,T) and P gas (n B ,T) are the energy den- 
sity and pressure of relativistic hadronic gas, respectively, 
which depend on baryon density n B and temperature T. 
The only difference from the ideal gas is that baryons are 
affected by a mean field U{n B ) = dW(n B )/dn B , i.e. the 
energy of the a-baryon of mass M a with momentum p 
is e a = (p 2 + A/, 2 ) 1 / 2 + b a U(n B ), where b a is the baryon 
number of the a-particle. The potential energy density 
W(n B ) is parametrized as follows 



W(n B ) 
= m N n 



+ 8,, 



a 



;as (n B ,0) 

\ 5/3 
n B \ 

n J 



\n J 



, (A3) 



where £ gas ("-s, 0) is the Fermi energy density of the cold 
nucleon gas 2 . This potential energy density depends only 
on the density n B . Parameters a, 6, and c are determined 
from the condition that the cold nuclear matter saturates 
at no = 0.15 fm -3 and s(no,T = 0)/na — mpf = —16 
MeV, and incompressibility of this nuclear matter is K. 
Basically we use the value K = 210 MeV (intermediate 
EoS). However, some calculations are also done for K — 
380 MeV (hard EoS), K = 130 MeV (soft EoS) and K = 
100 MeV (extrasoft EoS). 

Paramctrization i|A3(l results in supcrluminal sound ve- 
locity at high baryon densities. To preserve causality at 
high n B , the following form of the energy density 



s(n B ,T=0) = nompi 



n 



n B 



(A4) 



is used at n B /no > x c . Parameters A, B and C are 
determined on the condition that e(n B ,T = 0) and its 
two first derivatives are continuous at n B /no = x c . The 
baryon density, where paramctrization l)A4jl starts, de- 
pends on the incompressibility K: x c = 6 for our basic 
intermediate EoS, x c = 5 for the hard EoS, x c = 12 for 
the soft EoS. Parametrization (|A3(I for the extrasoft EoS 
does not violate causality at considered densities. 
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